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Q\ . ABSTRACT 
ON 

We discuss the present situation with regard to polarised nucleon structure 
Q H ' function measurements. In particular we examine the status of the Bjorken sum 

rule in the light of the recent data on the spin structure functions of (i) the 
deuteron obtained by the SMC group at CERN and (ii) the neutron by the E142 
group at SLAC. In order to fully exploit all the available data it is necessary to 
study the complete Bjorken system of equations, which may be done in any of sev- 
eral equivalent ways. Combining the new data with that already obtained for the 
' proton by the EMC group and earlier SLAC /YALE collaborations, together with 

£ — . bounds obtained on the strange quark polarisation, we show that the Bjorken 

system of equations is violated at about the 2a level. We also show, using un- 
polarised data, that arguments based on possible higher-twist contributions are 
unable to account for this discrepancy. In conclusion, a simple explanation in 
terms of a non-perturbative renormalisation of the relevant Wilson coefficients is 
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1. Introduction 

It has long been held (albeit by a limited community of "spin physicists") that 
polarised deep-inelastic scattering (DIS) can provide very stringent tests of factori- 
sation in perturbative QCD (PQCD) and the quark-parton model (QPM). Indeed, 
polarisation effects in general provide valuable insight into the dynamics of hadronic 
interactions and are extremely sensitive to the bound-state structure, so elusive to 
theoretical approaches.© In particular, the Bjorken sum rule (BSR)i'i is a measurable 
quantity that can be used for comparison between experimental data and theoretical 
predictions. The experimental precision attainable now surpasses the ten-percent level 
while, on the theoretical side, all relevant PQCD calculations have been carried out to 
at least two-loop order (i.e., one-percent level) and for the BSR itself to three loops.0 
Thus, one can consider such comparisons as serious, indeed obligatory, tests of the 
applicability of PQCD to such processes. 
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Like many other sum rules, the BSR is a direct consequence of the operator- 
product expansion, as applied to DIS and justified by asymptotic freedom in PQCD. 
In simple terms, all DIS structure-function sum rules may be derived from the hypoth- 
esized point-like behaviour of strong interactions at short distances. However, only the 
Adler sum rule! is directly associated with a conserved hadronic current and is there- 
fore unique in not receiving higher-order PQCD corrections, i.e., the Wilson coefficient 
associated with the relevant light-cone operator is unity to all orders in (and is in- 
dependent of) perturbation theory. That is, the validity of the Adler sum rule is not 
actually dependent on the validity of PQCD, being founded on symmetry principles of 
a deeper level, and its violation is virtually inconceivable in the present field-theoretic 
framework. In contrast, the Wilson coefficients associated with both the BSR and 
Gross-Llewellyn Smith sum rule! are 1— ^ (to first or next-to-leading logarithmic or- 
der); consequently, both sum rules are inextricably linked to perturbation theory and 
any discrepancy between their prediction and measurement should be taken as evidence 
for a failure of PQCD in predicting the normalisation of the Wilson coefficients. 

2. Polarised Deep-Inelastic Scattering 

By performing DIS with both beam and target polarised longitudinally, access is 
gained to the structure function g x {x, Q 2 ).i In the QPM this structure function has a 
simple relation to polarised quark distributions, analogous to that of Fi(x,Q 2 ): 

gi(x,Q 2 ) = ±E f e}Aq f (x,Q 2 ), 

F 1 (x,Q 2 ) = lE f e}q f (x,Q 2 ), 1 ' 

where the sum runs over quarks and antiquarks of flavour /, e/ is the /-quark fractional 
charge and x is the usual Bjorken scaling variable or momentum fraction of the struck 
quark. The quark densities are defined in the following manner: 

Aq f (x, Q 2 ) = qj(x, Q 2 ) - qj(x, Q 2 ), 
q f (x,Q 2 ) = qj(x,Q 2 ) + qj(x,Q 2 ), 

where qf(x, Q 2 ) are the densities of quarks of flavour / and positive or negative helicity 
with respect to that of the parent hadron. 

Of course, experimentally it is an asymmetry that is actually measured: namely, 
the ratio of the difference and sum of cross-sections for opposite helicity configura- 
tions. The polarised structure function is then extracted, using the following working 
definition: 

2 A(a;,Q 2 )F 2 (x,Q 2 ) 
^Q )= 2x(l + R(x,Q 2 )) ' (3) 

where Ri(x, Q 2 ) is the usual ratio of longitudinal to transverse unpolarised structure 
functions and Ai(x, Q 2 ) is the measured asymmetry. 
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Sum rules equate integrals over i B of such structure functions to independently 
known quantities. For example, using SU(2) current algebra and assuming scale in- 
variance, Bjorken showedlli that the proton-neutron difference for g x was given by the 
axial- vector /3-decay constant of the neutron, qa'- 

n~ n = f Q dxgr n {x,Q 2 ) = \g A . (4) 

Note that here there is apparently no longer any dependence on the energy scale, Q 2 . 
In fact, in PQCD there are radiative corrections to the right-hand side that depend 
on the running coupling constant as{Q 2 ) and thus only a logarithmic variation of the 
(experimental) left-hand side is expected at most. This sum rule therefore becomes 
a rigorous prediction of PQCD, through the justification of asymptotic freedom that 
leads to an approximate scaling behaviour. 

3. The Bjorken System of Equations 

The BSR, first dismissed as worthless! but later revalued,! is nevertheless difficult 
to test; it requires a precision measurement of g\ for both the proton and neutron over 
a sizeable range of xb (roughly speaking, a coverage of 0.01<a;<0.6 is necessary). While 
precise data have been available for the former target since 1988,01 the first ever data 
for the latter have only recently been published.&i!3 Moreover, the full SU(3) algebra 
of the baryon octet actually admits three independent quantities, which, while having 
their natural expression in terms of up, down and strange quarks, are better expressed 
in terms of the SU(3) axial-vector couplings: 

( I - ^7375^ I P T ) — 9Ai 

(P T I M7375M + ^7375^ - 2S7375S I p T ) = g A , (5) 
(P T I «7375« + d^y 3 j 5 d + S7375S \p 1 ) = 9o- 

The right-hand sides of the first two equations above correspond to known constants 
(c^=1.2573±0.0028@ and # A =O.629±O.O390), but the third (g ), corresponding to the 
flavour-singlet axial- vector current, is unknown. Thus an immediate prediction for, say, 
just the proton integral is not possible. 

Let us remark in passing that a further combination of the u, d and s axial- 
current matrix elements is, in fact, accessible in u-p elastic scattering^ and this would 
in principle allow an exact prediction for single nucleon targets. Unfortunately, the 
precision of such measurements is still too poor to permit any strong statement. 

Good arguments can be made, however, for setting the strange-quark matrix 
element equal to zeroS: there are very few strange quarks in the proton and they 
are concentrated below xb—0.1, where all correlations are expected to have died out. 
Therefore the second two matrix elements of eq. (^) were expected to be equal, thus 
leaving only two independent quantities. This allows predictions to be made for the 
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proton and neutron separately, which can now be expressed in terms of the two axial- 
vector constants and the strange-quark matrix element as 

rf n) = H^9A + h~9A + W I S "7375S I P r ), (6) 

For clarity, the PQCD corrections have been suppressed in the above, for more complete 
expressions see, e.g., ref. [L7|. Conversely, these equations may be used to extract the 
value of the strange-quark matrix element given the value of Ti for either nucleon, or 
(as in the case of the SMC) the deuteron. 

4. The Data — New and Old 

We shall now compare the experimental results obtained by the three experiments 
with the theoretical predictions based on the Ellis- Jaffe sum ruleS 
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-0.016 ± 0.006 ±| J As 



(8) 



The short-fall in the EMC proton measurement with respect to the Ellis- Jaffe prediction 
(taking /As=0) is immediately obvious. This observation led to the coining of the 
phrase Spin Crisis. A similar (though less striking) observation may be made for the 
SMC deuteron integral. In contrast, the neutron sum rule appears well satisfied by 
the E142 data. In terms of the strange-quark contribution, both the EMC and SMC 
measurements imply j'As~— 0.15 while that of E142 leads to JAs~— 0.02. Thus an 
important question is immediately raised: "How big can / As be?" This is the subject 
of the next section. 

Before moving on, let us make a few remarks on the SMC data: unfortunately, 
the errors are so large that they have little impact on any analysis although, taken at 
face value, they appear to agree with the BSR when combined with the EMC data. 
However, a certain neglect of these data can be motivated as follows. The particularly 
low value of the SMC deuteron integral is due, in roughly equal measure, to the large 
negative values reported at low x and to the very early drop at large values of x. 
While the low-x data are (barely) compatible with the smooth extrapolation to zero of 
the E142 data, the converse is certainly not true. Nevertheless, attempts of doubtful 
validity at combining the two data sets data have been made; the result is agreement 
with the BSR and (surprisingly) larger errors than the E142 data alone. This last fact 
clearly underlines the lack of statistical correctness of such a procedure. With regard 
to the high-x region, we have shown0 that the SMC data violate a rigorous bound 
that may be constructed from the combined EMC and unpolarised ,n data; the spin 
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data is used to eliminate the w-quark spin density and the F% ' n data then bounds that 
of the remaining rf-quark via positivity. Much of this discussion will be clarified by the 
planned E142 neutron measurements down to x~0.01, using a higher beam energy. 

5. Bounding the Strange-Quark Spin 

Positivity trivially implies \JAs(x)\<s(x), although we stress that equality would 
imply 100% polarisation. Now, Regge phenomenology tells us that the dominant con- 
tribution of the sea is is that of the pomeron (behaving as i" 1 ), which we recall is spin 
independent. We have used the very precise data on the unpolarised strange-quark dis- 
tributionsS to bound the non-diffractive contribution and thus too the strange-quark 
polarisation.Klro The result of this analysis is the following bound: 



As 



< 0.021 ±0.001, (9) 



which has been challenged (e.g., see refs. ^|,^5]) with various proposed constructions 
and/or parametrisations. However, we have shown in the papers cited above that all 
such proposals fail to agree with present experimental and phenomenological knowl- 
edge. 

A quantitative measure of the discrepancy between the data and theory may be 
obtained by performing a one-parameter fit for the strange-quark polarisation to the 
data, with and without the bound, the results^! are summarized in table [I]. 



Table 1. Results of a one-parameter fit for the strange-quark polarisation to the data, with 
and without the bound. 



Expt. Info. 


(p T I S7375S |p T ) 


X 2 /DoF 


C.L. 


EMC+Bound 


-0.027 ±0.019 


6.8/1 


0.9% 


EMC+SMC+E142+Bound 


-0.032 ±0.017 


10.6/3 


1.4% 


EMC+SMC+E142 


-0.088 ±0.020 


6.3/2 


4.3% 



6. Higher Twist and Higher-Order PQCD 

Not to be forgotten, of course, are the corrections to the various sum rules; as 
already mentioned above, there are PQCD higher-order corrections and the possibility 
of higher-twist contributions (especially in the case of the E142 data) should also be 
considered. For the non-singlet currents the PQCD corrections are known up to order 

However, since the value of as extracted experimentally is only valid up to order 
a s we only perform the analysis to this order. Although the inclusion of the second- 
order contribution does indeed shift the prediction for the BSR closer to the measured 
value, the effect is only of the order of a few percent. 
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The situation with regard to higher-twist contributions requires a little more care. 
Ellis and Karlinei@ have claimed that according to QCD sum-rule calculations^ the 
higher-twist contributions to the BSR actually bring the prediction perfectly into line 
with the measured value. In a similar analysis, Close and Roberts^ show that extremely 
large higher-twist contributions are required to explain all the data simultaneously. The 
argument is as follows: the neutron data are taken at an average Q 2 of only 2 GeV 2 
while the EMC data taken at an average Q 2 of 11 GeV 2 ; the problem then is to combine 
them at some common value of Q 2 . Experimentally it has been noted by each of the 
groups that the asymmetry is independent of Q 2 , within experimental accuracy. Thus it 
is reasonable to use the very precise parametrisations of F 2 (x, Q 2 ) to extract gi(x, Q 2 ) 
for any desired energy scale. Ellis and Karliner then propose to take the EMC data and 
"evolve" them down to the scale of the E142 data; it turns out that the variation in 
the integral of g\ is very small. At this point it is necessary to introduce a theoretical 
input, in the form of the above-mentioned predictions for the higher- twist component 
of the BSR. Such predictions are clearly model dependent and the associated errors are 
typically of the order of 100%. Moreover, a recent papeiHl has brought to light an error 
in the input for these calculations, suggesting that the contribution should be rather 
smaller and of the opposite sign. 

In any case, it is clearly preferable to try to avoid the the necessity of higher- 
twist estimates; this can be achieved by simply evolving the E142 data to the higher 
scales of the EMC experiment. It turns out, using this procedure, that the variation 
in the integral of is also small, reflecting a negligible higher- twist component. Thus 
the value of the BSR at Q 2 =\\ GeV 2 is only slightly affected by this procedure; the 
integrated values for the proton, neutron and BSR at Q 2 =\ \ GeV 2 are 

EMC r? = 0.128 ±0.019 

E142 Y\ = -0.032 ±0.011 (10) 

BSR (expt.) r?- n = 0.160±0.022 

to be compared with 

BSR (th.) rr n = 0.192±0.003 (11) 

Our results for the Q 2 dependence of the various sum rules are displayed graphically 
in figs. |1] and |2|. In table |2] we show the result of separating the various leading- and 
higher-twist contributions to the sum rules via a knowledge of the Q 2 dependence of 
the unpolarised structure functions. 

7. Extracting Spin Densities from the Data 

Adopting the attitude that the overall normalisation of the proton matrix ele- 
ments of the light-cone current operators is unreliably estimated in PQCD, it might 
be asked if there is any way to extract the three unknowns (i.e., Am, Ad and As) from 
essentially only two independent data sets (proton and neutron). Now, while overall 
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2 4 6 a in i ? 



Q 2 (GeV 2 ) 




Fig. 1. The solid curves are the Bjorken sum rule: a) asymptotic value, b) with the PQCD 
correction, c) PQCD and higher twist, d) experimental evaluation; the errors for a) and c) 
are given by the dashed lines. 



normalisations are often unreliable, ratios are well determined in the standard broken 
SU(6) picture of the nucleon; thus let us take as a reliable quantity only the ratio 

Au v ~ 2F ~ 3' 1 ' 

where F and D are the usual baryon-octet /?-decay constants; in SU(6) the ratio -§=§, 
experimentally it is ~O.6.0 Using this relation to substitute for Ad v in the expressions 
of for T"^ n and eliminating Au v between the two, we arrive at the following relation: 

r« = -^r? + §A s . (is) 

In obtaining the above, we have also included the phenomenological factor-two sup- 
pression of the strange over non-strange sea (although the final numerical results are 
rather insensitive to the actual value used). Thus, from our bound on As and the EMC 
value for Tf, we could have predicted 

0.002 > > -0.026, (14) 
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Fig. 2. The experimentally evaluated sum rules, as functions of Q 2 , of a) Bjorken, b) Jq 1 g\ 
and c) — J (solid curves) together with the PQCD (dashed curves) and PQCD plus 
higher- twist predictions (dotted curves). 

in rather good agreement with the SLAC data. 

Alternatively, we can turn these expressions around and extract the various inte- 
grated quark spin-dependent densities, using both EMC and SLAC data: 

Am, = 0.72 Au = 0.66 

Ad v = -0.24 Ad = -0.30 

Au s = Ad s = -0.06 As = -0.03 ( 15 ) 

AS = Au+Ad+As = 0.33 

(the errors on the above quantities are typically of the order of 0.02 — 0.03). Notice 
(i) the strange-quark component is entirely compatible with the bound and (ii) the 
total quark spin, AE, is a sizeable fraction of that of the proton. With the additional 
standard assumption of an extra factor of 1—x for the valence d density with respect 
to u, the distributions as functions of x may be obtained; the corresponding curves are 
shown in fig. |^. 
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Table 2. The proton and neutron integrals and higher- twist corrections, the fourth and fifth 
columns give the two corrections and the last column contains the corrected (i.e., leading- 
twist) integrals. Note that rounding errors cause slight discrepancies between the various 
figures. 
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-0.001 


0.000 
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0.128 


0.000 


0.000 


0.127 ± 0.017 


p — n 


0.159 


0.000 


0.000 


0.159 ±0.020 



8. Conclusions 

Let us preface the closing remarks by noting that all of the above is crucially 
dependent on the validity of the original EMC data, which we have taken at face value. 
Thus, given the enormity of the conclusions to which we are ineluctably drawn, it is 
important to stress the urgency of remeasuring the proton spin sum rule to the highest 
possible precision. 

Although the BSR only fails to be satisfied experimentally at something like a \\a 
level, the complete system of equations including the strange-quark polarisation bound 
is violated at a level of 2|er. We have also shown that a mild assumption of incorrect 
evaluation of the overall normalisation for hadronic matrix elements in PQCD leads to 
a picture of the nucleon in which a large fraction of the spin is carried by the quarks 
and the strange-quark spin is very small, i.e., precisely the picture existing before the 
beginning of the Spin Crisis. 

To conclude, let us just recall that a correct prediction of the proton spin struc- 
ture function had appeared in the literature^ in 1985, well before the EMC results. 
Implicit in the ACD Fire-String model that led to this prediction, is the almost vanish- 
ing neutron spin asymmetry, as corroborated by the E142 results. The model, with an 
asymmetry driven essentially by a broken SU(6) symmetry, has a final-state structure 
in which non-physical (i.e., coloured) states are explicitly barred from contributing. 
The two contributions to DIS are (i) the single fire-string diagram, where quark he- 
licity conservation makes itself strongly felt and the broken SU(6) structure leads to 
essentially vanishing neutron asymmetry; and (ii) the double fire-string diagram where 
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Fig. 3. Our extraction of the quark spin distributions: the upper and lower dashed lines are 
the u and d quarks respectively, the dotted line is the s quark and the solid line is the resulting 
total quark contribution. 

the effects of quark polarisation are so dilute as to be negligible in all cases. The differ- 
ences between the explicitly physical final-state structure of ACD and that implicitly 
unphysical in PQCD explain the failure of the latter to obtain the correct normalisation 
(or Wilson coefficients) of the nucleon matrix elements in question. 

9. Note added 

Since presenting this talk the next-to-leading PQCD corrections to the Ellis-Jaffe 
sum rule (in other words for the flavour- singlet piece) have been made available.il 
Thus, given the value of any single nucleon integral at some Q 2 , one can perform the 
two-loop evolution to any desired Q 2 scale. We have checked the effect of this approach 
and find negligible difference with the results presented above, although it should be 
remarked that the theory-data discrepancy is then, in fact, slightly accentuated. 
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